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J^ ■ Abstract 



We study spatial correlation functions of charmonium in 2+1 flavor QCD using an improved 
staggered formulation. Contrary to the temporal correlation functions the spatial correlation func- 
tions exhibit a strong temperature dependence above the QCD transition temperature. Above this 
_^ temperature they are sensitive to temporal boundary conditions. Both features become significant 

(— I . at a temperature close to 1.5Tc and suggest corresponding modifications of charmonium spectral 

functions. 
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I. INTRODUCTION 



Quarkonium suppression was proposed long ago by Matsui and Satz as a signal for decon- 
finement in heavy ion collisions [l| . The basic idea behind this proposal was the fact that at 
sufficiently high temperatures color screening will prevent formation of heavy quark bound 
states (for recent reviews on this topic see [2, |3[). Charmonium suppression was indeed 
observed in heavy ion experiments at SPS ^ and at RHIC [5|, and is now being intensively 
studied also at the LHC 6|]. In order to interpret current and future experimental findings 
it is very important to know (among other things) the in-medium properties of heavy quark 
anti-quark pairs. 

The notion of bound state melting can be rigorously formulated in terms of the spectral 
functions. Charmonium dissociation corresponds to gradual broadening and eventual disap- 
pearance of the bound state peaks in the corresponding meson spectral functions. Spectral 
functions are related to the Euclidean time meson correlation functions and can be studied 
using lattice QCD calculations. The standard approach to obtain information about the 
spectral functions from calculations of temporal correlators relies on the Maximum Entropy 
Method (MEM) 7|, |8|. Based on this approach, analyses in quenched QCD led to the con- 
clusion that IS charmonium states may survive up to temperatures as high as 1.6Tc, with Tc 
being the deconfinement temperature of the SU(3) gauge theory 9|-[l2| (for a similar analysis 



in 2 fiavor QCD see 13|). Other lattice QCD studies based on different techniques, e.g. vari- 
ational analysis, led to similar conclusions IJ, ll5[ . The conclusion that charmonium states 
survive in a certain temperature range above the QCD transition temperature is largely 
based on the observation that the corresponding temporal meson correlation functions show 
a weak temperature dependence across the transition jlll . 
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16 



4l8l|. Moreover, at high 



temperatures in particular the vector spectral function receives a contribution near zero 
frequency, a; ^ 0, which corresponds to the transport of heavy quarks. This near-zero-mode 
of the spectral function gives rise to a contribution to the temporal correlation function 
that is (almost) constant in Euclidean time and is responsible for most of the temperature 



dependence of temporal meson correlators [l6|, ll7|. Thus, the temperature dependent effects 



due to in-medium modifications and dissolution of charmonium states appear to be small in 



temporal correlation functions ^ 

Contrary to the lattice analysis of temporal meson correlators, potential model studies 



that use the static quark anti-quark correlators calculated in lattice QCD 20| as an input 



into the Schrodinger equation, predict melting o 
slightly above the QCD transition temperature 



charmonium bound states at temperatures 



21 



23| . Furthermore, analyses of quarkonia 



at non-zero temperature within the effective field theory approach revealed that in addition 



to the modification of its real part, the potential also acquires an imaginary part [24, |25 |. 
The imaginary part plays an important role in quarkonium dissociation. Even when one 
takes into account the uncertainties in relating the imaginary part of the potential to the 
correlation functions of a static quark anti-quark pair calculated in lattice QCD, a non- 
vanishing imaginary part leads to the dissolution of at least the charmonium bound states 



26|. A similar conclusion was reached in the analysis of the charmonium spectral functions 
using the T- matrix approach [27|. Remarkably though, it was shown in the framework of 
potential models that the melting of quarkonium states does not result in large changes of 
the temporal correlators and the correlators seemed to be temperature independent to very 



good approximation [22|. Thus both potential models and direct lattice QCD calculations 
suggest only a weak temperature dependence of the temporal meson correlators across the 
transition. 

In order to extract detailed features of the spectral functions, e.g. broadening of bound 
state peaks or melting of states, from such weakly temperature dependent temporal corre- 
lators one needs high precision data at a large number of Euclidean time separations. The 
most recent high statistics (quenched) lattice calculation of charmonium spectral functions 



using MEM reached lattice spacings down to a = 0.01 fm [18|. This study found no evi- 
dence for bound state peaks at T ^ 1.46Tc, suggesting charmonium melting to take place at 
a temperature somewhat smaller compared to the earlier MEM based lattice QCD analy- 
ses. Nevertheless, it is desirable to have alternative observables that may provide additional 
information on the in-medium properties of heavy quark bound states. 

In this paper we therefore analyze spatial charmonium correlators at non-zero tempera- 



ture. Following 



we also study the dependence of the spatial meson correlation functions 



^ Using NRQCD at finite temperature some evidence for melting of bottomonium P-states, that have same 
size as the ground state charmonium, was obtained [l9|. 



on the temporal boundary condition. We have found that the behavior of spatial char- 
monium correlators and their dependence on the boundary condition is consistent with 



significant modifications of charmonium states at sufficient 



liminary results of this work were presented previously in 29 



1. 



ligh temperatures. Some pre- 



30| . We also note that some 



exploratory studies of spatial charmonium correlators in quenched QCD were reported in 



ll| . Our findings in 2+1 fiavor QCD are qualitatively consistent with these calculations. 



II. SPATIAL CORRELATION FUNCTIONS 

Correlation functions of a meson operator J = qTq along the spatial direction z are 
defined as 

G{z, T)= dxdy I dT{J{x, y, z, r) J(0, 0, 0, 0)). (1) 

The matrix F is a product of Dirac matrices and fixes the quantum number of the meson. 
The spatial correlation function is related to the meson spectral function at non-zero spatial 
momentum 

G{z,T)= / dp,e^P^V(a;,p„T) (2) 

Jo ^ J~oo 

Thus the temperature dependence of the spatial correlation function also provides infor- 
mation about the temperature dependence of the spectral function. Medium effects are 
expected to be the largest at distances that are larger than 1/T. At these distances G{z,T) 
decays exponentially and this exponential decay is governed by a screening mass Mgcr- It is 
also pertinent to note here that the above relation suggests, unlike the temporal correlation 
functions, an uj6{u) type zero- mode contribution to the spectral function does not lead to a 
non-decaying constant contribution to the spatial correlation functions and only leads to a 
contact term. 

If there is a lowest lying meson state of mass M, this is signaled by a peak in the spectral 
function, a{u,Pz,T) ~ 6{uj'^ ~ vI ~ M"^). This bound state peak in the spectral function 
determines the long distance behavior of the spatial meson correlation function. We thus 
have Mgcr = M. On the other hand, at very high temperatures the charm quark and anti- 



quark are not bound and the meson screening mass is given by 2 A/(7rTp + ?7if , where rric 
is the charm quark mass and vrT is the lowest fermionic Matsubara frequency. The above 
value of the screening mass is a consequence of the anti-periodicity of the fermion fields 
along the time direction and forces the charm quark (anti-quark) to pick up at least a ttT 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



G(zJ)/G(z,T=0) 


















■ 1 - - . 

• ♦ 


• 
* 


• 


• 


• 


• • 


• 


i 


* i i 


• 




• 


• 












• 








• 












• 


• 






• 
• 


• 


• 










• 










• 


T/T(, =1 .23 






H 


• 


^ 






• • 




"■ 


1.84 
2.44 




















* 


H 




• 


* 


# 


- 






















• . 


1 








z 


[fm] 






• • 



0.0 0.2 0.4 0.6 



0.8 



1.0 1.2 



FIG. 1: Ratio of the spatial pseudo-scalar correlators to the corresponding zero temperature cor- 
relators calculated for a~^ = 2.8 GeV at three different temperatures. Calculations have been 
performed on lattices with temporal extent A'^^ = 6, 8, and 12. In units of the transition tempera- 
ture, these correspond to T/Tc = 1.23, 1.84, and 2.44. 



contribution from the lowest non-vanishing Matsubara frequency at non-zero temperatures. 
The transition between these two hmiting behaviors of the screening masses may serve as 
an indicator for significant modifications and ultimately dissolution of the meson states. 

Furthermore, if the cc pair forms a mesonic bound state the screening mass is not expected 
to be sensitive to the anti-periodic temporal boundary condition due to the bosonic nature 
of the basic degrees of freedom. In this case using a periodic temporal boundary condition 
for the fermions is expected to give the same screening mass as for the case of an anti- 
periodic temporal boundary condition. Since for the non-interacting theory with a periodic 
temporal boundary condition Mgcr = '^f^ic, a comparison with the usual anti-periodic case 
of Mscr = 2y(7rT)2^rm| will also facilitate in the identification of free theory like behavior 
of the spatial correlation functions. Thus studying the dependence of the screening masses 
on the temporal boundary conditions may also provide some additional information about 
the existence of bound states. 



III. NUMERICAL RESULTS 



We performed numerical calculations in 2+1 flavor QCD using the improved staggered 



p4 action 32 



33| on 32'^ x A^^ lattices with A'^^ = 6, 8, 12 and 32. Some calculations have also 
been performed on 24^ x 6 lattice. We mostly used the gauge fleld conflgurations generated 



for the study of the QCD equation of state 3J, |35|. The strange quark mass rris was flxed 



to its physical value, while for the light quark masses we used rris/lO corresponding in the 
continuum limit to a pion mass of about 220 MeV. The lattice spacing was flxed using the 



value vq = 0.469 fm 



the lattice parameters is given in Refs. 



M 



for the Sommer scale [37[. A detailed discussion of the choice of 
35|. Since for lattice spacings used in this study 



cutoff effects may be still signiflcant we will present our results in terms of the reduced 



temperature T/Tc. For the transition temperature 
values T, = 198 MeV and T, = 191 MeV from Ref. 



or Nt = 6 and 8 lattices we will use the 



38 



39| . From the 0{N) scaling analysis 



presented in 39| we estimate Tc ~ 160 MeV in the continuum limit for light quark masses of 
rris/lO. Assuming a 1/N^ dependence of T^ for A^,- > 8 we estimate that T^ ^ 174 MeV for 
Nt- = 12. We also used the staggered p4 action for the valence charm sector. The valence 
charm quark masses were flxed to reproduce the physical value of the J/ip mass at T = 



4G 



4l|. 



The staggered fermion formulation describes four valence quark flavors in the continuum 
limit. Meson operators in this formalism are written as J = q(T^ x T^)q, where F^ and F^ 
are products of the Dirac matrices and describe the spin flavor structure of the meson 42|. 
Here we consider only local meson operators for which F^ = F^ = F, since these are the 
cheapest to calculate. As in most lattice QCD studies only the quark-line connected part of 
the correlation function is taken into account; however, for heavy quarks the effect due to 
the disconnected part is expected to be small. The meson operators can be written in terms 
of staggered quark flelds, x ^^^ x ^s J = (j){X)x{X)x{X), X = {x, y, z, r) |43|, |4^. Different 
meson channels/operators are flxed by appropriate choice of the phase factor </)(X). Each 
choice of 4>{X) corresponds to a pair of F matrices related to positive and negative parity 
states that show up as oscillatory and non-oscillatory terms in the correlation function 44 1. 
In our analysis we used both point and wall sources. 

In this paper we mainly discuss the simplest case corresponding to (piX) = — (— 1)"^+^+^. 
We will comment on other choices of (f){X) later. The correlator corresponding to (pi^X) = 



6 



— (—1)'^"'"^"'"^ is dominated only by pseudo-scalar meson states (the oscillatory contribution 
vanishes) and thus is the easiest to analyze. For this reason we refer to the corresponding 
correlator as the pseudo-scalar correlator. In the zero temperature limit this correlator 
corresponds to the spin singlet S-wave charmonium state rjc- First we are interested in the 
change of the correlation functions in the high temperature regime of QCD relative to the 
zero temperature correlator. In Fig. [T]we show the ratio of the spatial correlators calculated 
for three different temperatures corresponding to A't- = 6, 8 and 12 and for fixed lattice 
spacing a~^ = 2.8 GeV. As one can see from the figure the correlator shows significant 
change already at a temperature of T/T^ ~ 1.2. As the temperature increases further the 
deviations from the zero temperature correlator become more prominent. This is different 
from the case of the temporal correlators, where very little change is seen even at the 
highest temperature. This is partly due to the fact that larger separations can be probed in 
the spatial direction. Indeed, for separations smaller than 1/(2T), available in the case of 
temporal correlators, the temperature dependence is very moderate. 

We have studied the long distance behavior of the spatial correlators and using simple 
exponential fits extracted the corresponding screening masses as well as the amplitudes. In 
the zero temperature limit the latter are proportional to the square of the wave functions 



at the origin. At infinite temperature, i.e. in the 



ree case limit, the leading behavior of the 



correlation function has the form ~ e~"^'"'''^/z [45|. However, a 1/z dependence of the am- 
plitude or, corrolary, a decrease of the effective screening mass rriscriz) = \nG{z)/G{z + 1) 
proportional to ln(l -|- 1/z) ~ 1/z has not been observed. 

In Fig. [2] we show the screening masses and the amplitudes divided by the corresponding 
zero temperature values as function of temperature in units of the transition temperature 
Tc. As one can see from the figure the cutoff effects for the screening masses and the 
corresponding amplitude are small. At low temperatures both the masses and the amplitudes 
agree with their T = values within errors. Above the transition temperature we see small 
but statistically significant deviations from the zero temperature result. We also compare 
our numerical results for the screening masses for T > 1.2Tc with the free theory result 



2y(7rTp + mf. Here we used m^ = 1-42 GeV, as will be justified later. 

For T > 1.5Tc deviations from the zero temperature limit become large and the screening 
masses are compatible with the behavior characteristic for unbound cc pairs. The behavior 
of the charmonium screening masses is very similar to the findings of a study in quenched 
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FIG. 2: Pseudo-scalar screening mass (a) and amplitude (b) divided by the corresponding zero 
temperature values as function of the temperature. Filled symbols show the results obtained with 
point sources, while open symbols represent the results from the wall sources. The quenched results 
111 ] and correspond to the smallest lattice spacing [pT^ '^ 9.7 GeV) used in 



are taken from Ref. 



that study. (This corresponded to lattices with temporal extent N^- = 12, 16, and 24). The solid 
black line is the free quark result (see text for details). 



QCD ll|. The change in the temperature dependence of the charmonium screening masses 
around T = 1.5Tc is similar to the change of the screening masses in the light quark sector 
around the QCD transition temperature 3l|- Therefore the change in the behavior of the 
charmonium screening masses around T = 1.5Tc is likely due to the melting of the meson 
states. 

As stated above, the dependence of the screening masses on the temporal boundary 
conditions can provide additional information about the existence of bound states in the 
deconfined medium. Therefore, we also calculated the pseudo-scalar charmonium screening 
masses using periodic boundary conditions in the time direction. The comparison of the 
screening masses calculated with periodic and anti-periodic boundary conditions is shown 
in Fig. 131 As one can see from the figure there is no dependence on the boundary conditions 
in the low temperature phase of QCD. This is consistent with the fact that cc pairs form 
bound states. For T^. <T < 1.5Tc we start seeing small differences in the screening masses 
calculated with periodic and anti-periodic boundary conditions. This might be an indication 
of broadening of the IS charmonium state. At temperatures T > l.ST^ we already see a 
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FIG. 3: Pseudo-scalar charmonium screening masses calculated with anti-periodic (filled symbols) 
and periodic (open symbols) boundary conditions as function of the temperature. Solid (dashed) 
lines correspond to the free theory prediction with nic = 1.42 GeV for anti-periodic (periodic) 
boundary conditions. 

strong dependence on the boundary conditions indicating that the fermionic sub-structure 
of cc pairs becomes relevant and the charm quark and anti-quark start becoming sensitive 
to the lowest Matsubara mode, which is non-zero in the case of anti-periodic boundary 
conditions. This is expected to happen when the bound state is dissolved. If periodic 
boundary conditions are used charmonium screening masses are equal to 2mc in the limit of 
very high temperature. If we choose rric = 1.42 GeV at the highest temperature we recover 
this expectation. Therefore, we used the value rric = 1.42 GeV when comparing with the 
free field theory prediction. Note, however, that for T > 2Tc the dependence of the screening 
masses on the boundary condition is compatible with the free theory expectation regardless 
what value of rric is chosen; i.e. we can convert between the results obtained with periodic 
boundary condition to the ones obtained with anti-periodic boundary condition by simply 
adding the lowest Matsubara frequency in quadrature. 

Let us comment on the behavior of the screening masses in other channels, i.e. for other 
choices of 0(X). The choice 0(X) = — (— l)^"*"* (or 0(X) = — (— 1)^"*"*) corresponds to F = 71 
(or F = 72) in the non-oscillatory and F = 7270 (or F = 7170) in the oscillatory part, 
which in turn correspond to J/ip and he, respectively. We have found that the correlator 
for this meson channel is dominated by the J/ip and we were unable to obtain any signal 



9 



for the he- The quahtative 



similar to that of the rjc 29 



Dehavior of the J/ip screening masses was found to be quite 
30|. Finally, for F = 1 a pseudo-scalar {r]c) and a scalar (xco) 
states contribute to the non-oscillatory and oscillatory parts of the correlator. Similarly the 
choice (f){X) = —(—1)'^ (or (j){X) = —(—I)''' ) corresponds to 7273 (or 7173) and 7275 (or 
7175). Thus vector {J/ip) and axial- vector (Xci) states correspond to the non-oscillatory 
and oscillatory parts of the correlator. Contrary to the situation with light quarks, all three 
channels remain to receive both, non-oscillatory and oscillatory contributions above Tc. In 
fact, up to temperatures of about 2Tc the non-oscillatory part is larger. Thus the screening 
masses corresponding to the non-oscillatory part of the correlators can be reliably calculated 
and show a temperature dependence that is very similar to that of the pseudo-scalar (77^) 
and vector (J/ip) screening masses discussed above. On the other hand, the scalar and 
axial-vector screening masses could be reliably determined only at temperatures T > 2Tc. 
At these temperatures they are consistent with the free theory value. 

IV. CONCLUSIONS AND OUTLOOK 

In this paper we have demonstrated that spatial charmonium correlators are useful tools 
for studying in-medium charmonium properties and their possible dissolution in the decon- 
fined phase. The pseudo-scalar screening masses become temperature dependent just above 
the crossover temperature of 2 + 1 flavor QCD indicating that medium modifications of the 
charmonium spectral function set in gradually. At temperatures T > 1.5Tc we see already 
large modifications of the charmonium correlators, and the temperature dependence of the 
screening masses is qualitatively the same as the behavior expected for unbound cc pairs. 
At T > 2Tc the value of the screening mass and its dependence on the boundary conditions 
is in agreement with the free field limit indicating the absence of charmonium states at these 
temperatures. The lattice results presented here are consistent with the expectation based 
on potential model and can be used to further constrain potential model calculations of 
the charmonium spectral functions. This, of course, will require that those calculations are 
extended to finite spatial momentum. 

The investigations presented in this paper can be extended in different ways. While 
above the transition temperature cutoff effects are under control, this is not the case in 
the hadronic phase. To reduce the cutoff effects at low temperatures calculations using the 
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Highly Improved Staggered Quark (HISQ) 46|| action will be necessary. The HISQ action 



can also reduce cutoff effects in the light quark sector. It would be furthermore interesting 
;o do the calculations with a Wilson type action for heavy quarks, like the Fermilab action 
47|. Since in this case scalar and axial- vector do not have to be disentangled from S-wave 
states this would presumably allow for a more accurate study of the P-wave charmonia. 
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